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SUMMARY 

Active magnetic bearings used in rotating machinery should be designed as lo- 
cally controlled, independent devices similar to other types of bearings. The func- 
tions of control electronics and power amplifiers can be simply and explicitly 
related to general bearing properties such as load capacity, stiffness, and damping. 
In this paper, the dynamics of a rotor and its supporting active magnetic bearings 
are analyzed in a modified conventional method with an extended state vector con- 
taining the bearing state variables. 


INTRODUCTION 

Active Magnetic Bearings (AMBs) have been slow to show acceptance in the rotat- 
ing machinery industry, and to machinery designers, they remain somewhat a mystery. 
AMBs have been developed over the decades by researchers in electrical engineering 
and systems control. Consequently, progress in their development is often presented 
in a language that is alien to the mechanical designers. The design engineer, whose 
major discipline is mechanical engineering, is likely to be unfamiliar with the con- 
cept that AMB stiffness and damping are functions of excitation frequency. In addi- 
tion, the designer does not commonly use terminology such as control bandwidth, 
phase compensation, control spillover, etc. Not only does the AMB's electrical pre- 
sentation hinder its acceptance, but also the fact that it is such a small component 
of a sophisticated machine; most design engineers will question whether it warrants 
the time and effort needed to gain familiarity with it. Consequently, greater un- 
derstanding of the AMB is imperative in increasing its use In conventional 
rotor-bearing applications. With these observations in mind, the authors herein at- 
tempt to describe the key elements and properties of AMBs in a general engineering 
language that includes rotor dynamic notation which mechanical engineers are more 
accustomed to using. 


BASIC PRINCIPLES 


The principles of an AMB are well explained by Habermann and Liard (1980). The 
practical AMBs generally adopt an 8-pole stator configuration as shown in Figure 1. 
Both the stator and journal are made up of stacks of laminations of ferromagnetic 
material. The journal is shrunk on a shaft without windings. Laminations reduce 
eddy currents that not only create a power loss but also degrade the performance of 
the bearing. The eddy currents generate magnetic fields that oppose the flux needed 
for regulating the journal motion. The eight poles of the stator are separated into 
four quadrants (Weise, 1985). In each quadrant, the electromagnetic windings are 
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wound in such a way that the magnetic flux will circulate mainly inside the quad- 
rants (Chen and Darlow, 1987). In other words, the magnetic forces of the poles in 
any quadrant can be varied by changing the amount of current in the windings of that 
quadrant without causing force changes in other quadrants. Therefore, each quadrant 
of poles can be controlled independently. 

It is well known that the magnetic force is proportional to the current to air 
gap ratio squared (see Figure 2). Engineers conscientiously try to avoid nonlinear 
design and analysis. The linearization of AMB dynamics is achievable by making the 
air gap large relative to the journal normal excursion and by providing a relative- 
ly large steady state current (called bias current hereafter) through each quadrant. 
The bias current produces I^R loss, which is a major power loss in AMB. However, 
the total resistance including the windings and other parts in current path is not 
large; the AMB power loss in general is insignificant when compared to the conven- 
tional bearings, such as the hydrodynamic oil-film type. 

It is noteworthy that all poles exert attractive forces on the journal. Two 
opposite quadrants of poles are used to center the journal in one direction. While 
generating cancelling forces may seem to waste electric power, the AMB stiffness and 
damping are directly proportional to the bias current. 

The journal floating in a magnetic field that is produced only by the bias cur- 
rents is not stable. This situation is analogous to supporting a vertical stick at 
the bottom by a hinge without controlling the hinge. To create stability at the AMB 
center, the journal motion must be sensed and corrected instantaneously and contin- 
uously by superimposing a small control current to each bias current. For example, 
when the journal moves upward off the center by a small displacement, Y, the current 
in the top quadrant will be reduced by a small amount, i, and the bottom quadrant 
increased by i. The control currents produce a net downward force (-F) that pulls 
the journal back to the center. From sensing Y to producing -F, a series of AMB 
components become involved. It is the authors 1 preference to separate them into two 
groups : 


• Sensor and control electronics 

• Power amplifiers and electromagnets. 

Before examining them in detail, presenting basic statics and dynamics will 
specify the function of these components. 

AMB STATICS AND DYNAMICS 

To support a static load, W, in Y-direction, the net attractive force caused by 
the bias currents is 

W = fdx 2 - I 3 2 )/C 2 (1) 

where f is a magnetic pole constant for a given number of windings. 

Choosing an I 3 , I 3 can be determined by equation 1. I 3 must be larger than the 
control current, i, at any anticipated transient situation so that no current satu- 
ration can occur. 

Assuming that 

If, I 3 » i and Y « C 
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and that the AMB sees a dynamic mass, M, then a linearized equation of motion is 

M V - K m Y = Ki i (2) 

where K[ and K m are respectively called the current stiffness and the negative 
spring of the magnetic field due to the bias current* Without the small control 
current, i, equation 2 represents an unstable dynamic system. To make it stable, it 
is simple and sufficient to have 

i = G E 

and 

E = -C d Y - C v Y 

where E is the controller output voltage signal and the input to the amplifiers, and 
G is the power amplifier gain (amp/volt). In other words, the small control current 
should be made proportional to the journal displacement and velocity. The minus 
sign is indicative of negative feedback, which is a common regulating mechanism. C d 
and C v are constants, and called proportional and derivative control gains. 

Substituting equations 3 and 4 into 2 produces 

M Y + GKiC v Y + (CK^Cd - Km) Y = 0 (5) 

If C d is large enough, such that GK[C d > K m , we achieve a stable one degree of 
freedom system in Y-direction which has damping coef f icient GK£C V and stiffness coef- 
ficient (GK£C d - K m ). 

If the static load W, is miscalculated, or there is an occasional slowly vary- 
ing load, (as can be caused by the gyroscopic effect of a manuvering spacecraft) the 
journal will sag or drift off the bearing center. To correct this potential prob- 
lem, a third corrective mechanism called the integral control, as depicted in Figure 
3, can be added to equation 4. That is, 

E = -C d Y - C v Y -C e J Ydt (6) 

As shown in Figure 3, the accumulated journal position error over a period of 
time will produce a part of the control voltage. The integral control does not re- 
spond to high frequency vibrations. 

SENSOR AND CONTROL ELECTRONICS 

For AMB journal displacement measurements, three practical sensors exist: the 
capacitance probe, inductance probe, and eddy current probe. Each has advantages 
and disadvantages, but they all relate the small distance (in mils), from the sta- 
tionary sensor to the rotating shaft, to an output electrical signal in volts. 

No inexpensive and reliable sensors are available yet for measuring the journal 
velocity (Y). Consequently, different control circuit designs have been developed 
to produce a psuedo velocity from the displacement measurements. The feedback of 
velocity is in essence a corrective action of anticipation. If the vibration dis- 
placement is sinusoidal, the corresponding velocity is also sinusoidal, 

Y = Y 0 sin(wt) 

Y = U)Y 0 sin(ut + tt/2) 


(3) 

(4) 
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but peaks 90° earlier on an oscilloscope. The logical way of obtaining the velocity 
from the measured displacement is by performing differentiation with an analog cir- 
cuit. However, inevitable error or noise signals occur in this measurement; a small 
displacement error at high frequency will become a large erroneous velocity at the 
same frequency by differentiation, which may overheim or bury the real velocity sig- 
nals. Therefore, a practical differentiator always applies a certain M low-pass fil- 
tering 11 , which reduces the high frequency signal content. A conflict arises between 
reducing the high frequency signal and keeping the lead time for the low frequency 
signal. The filtering aspect of circuit design may become awkward at times. 

A more popular method of generating the velocity component of controllers is to 
use an analog circuit called phase-lead network. Each phase-lead network provides 
lead time or phase advance for the displacement signal in a reasonably large fre- 
quency range. The maximum phase advance per stage network is less than 70°. Two 
networks may be needed to cover a large spread of critical speeds. The output of a 
phase-lead network is a combination of displacement and velocity signals. Thus, its 
feedback produces the damping and contributes to the stiffness. 

An alternative method, called the Velocity Observer, was recently developed and 
used by Chen and Darlow (1987) for producing the appropriate control signals. It 
integrates journal force, (equivalently, acceleration) to obtain velocity. The ob- 
server takes the displacement and the small control current as input. The Velocity 
Observer has an advantage over the phase-lead approach in that it covers much wider 
frequency range with phase advance > 90°. 

The integral feedback term in equation 6 is straight forward in analog circuit 
design. The integral control is a time delay action. To confine its effect below 1 
or 10 Hz, where system resonance seldom occurs, low-pass filtering may be included 
in the integrator circuit. A comparison of the mathematical expressions of the 
three types of circuits are presented in Figure 4. 

The control electronics, which take the displacement Y as input, is generally 
called a Proportional-Integral-Derivative (PID) controller. There are many ways to 
implement the controller. Humphris et al. (1987) used a single path proportional 
derivative (PD) controller with two phase-advancing circuits in series. Fukata et 
al . (1986) used a PID controller with three parallel paths and a differentiator. 

Chen and Darlow (1987) used another PD controller using two parallel paths and a Ve- 
locity Observer. 

The parallel-path approach is preferable because the adjustments of the stiff- 
ness and damping are less dependent on each other. A controller example is shown in 
Figure 5 using a single stage of phase-lead network. The transfer function of the 
controller can be expressed in the following closed-form formula: 

-E/Y = C d + C v (s + a)/(s + b) + C e /(s + Wj) (7) 

From (2), the AMB dynamic force is 

F = K£ i + K m Y. (8) 

Assuming the amplifiers provide control current as much as demanded with no delay, 
and for simplicity G = 1 (amp/volt), then 

i * E (9) 
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Combining equations 7, 8, and 9, we have (for s - jw) 

-F/Y = K + jcd B (10) 

where ? t ? o ? 

K = Ki[(C d -K m /Ki)+C v (ab+co z )/(b z +U) z )+C e W 1 /(w 1 z +a)^)] (11) 

B = K i [C v (b-a)/(b 2 +w 2 )-C e /(a) 1 2 +w 2 )] (12) 

Equations (11) and (12) show that AMB stiffness and damping are functions of excita- 
tion frequency, U). An example of the stiffness and damping using this controller 
is presented in Figure 6. It may appear unusual that bearing damping in Figure 6 is 
negative in the low frequency range. The fact is, that as long as no system vi- 
bration mode exists in that range, there should be no instability. Note that for 
every U), a damped frequency value oj can be calculated from the bearing parameters 
M, K and B. A natural vibration mocTal frequency exists at W = 

POWER AMPLIFIERS AND ELECTROMAGNETS 

The physical size of an AMB journal and stator can be determined (Chen 1988) by 
first calculating the pole surface area: 

Ap = Fmax/250 (13) 

The formuLa was derived in American units for silicon steel laminations (satu- 
ration flux 110 x 10~* lines/in. 2 ), A p is taken as the smallest cross-sectional area 
(in, 2 ) along the pole. F max is the maximum force (lb) that two opposite pairs of 
poles can take without flux saturation. Choosing the axial length L p , the circum- 
ferential pole width is A p /L p . Then, the radial dimensions can be decided outward 
from a given shaft diameter at the AMB. The sizing guidelines are 

• The cross-sectional area at any point of the flux path is not less than A p . 

• Adequate wiring space and cooling surface are provided. 

• The axial length is no greater than the journal outer diameter. 

• As a rule, the air gap should be ten times the expected journal excursion. 

The smaller the gap, the less power is required. 

Converting a low power control voltage signal to a high power control current 
and actuating the electromagnets requires current source or power amplifiers. Two 
types of power amplifiers, the linear type and the pulse-width-modulation (PWM) 
type, have been studied and used by MTI in their magnetic bearings. The linear type 
applies the control signal to a power transistor in "active" mode. The transitor 
continuously regulates the current through the windings from a DC source. In this 
mode of operation, the voltage drops across the transitor, and thus the power loss 
is high. The PWM type applies the control signal to generate high voltage output 
pulses at a frequency above audible range. The on-time period of each pulse is pro- 
portional to the input signal. The high voltage pulse train produces current in the 
windings. The PWM amplifier is electrically noisy but much more efficient in pro- 
viding the required power with low power loss. The power transitors in it operate 
in "saturation" mode and have small resistances. Schweitzer and Traxler (1984) had 
indicated that the borderline in favoring one type over the other is about 0.5 kVA. 

The inductance of electromagnet windings can cause a problem of control current 
delay. The inductance, L, is proportional to N^A p /C. A time constant, L/R, is as- 
sociated with the power amplifier-electromagnet system. A controller signal at w > 
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R/L will have a phase delay > 45°, which is one-fifth of a cycle. Also, the induc- 
tance limits the permissable rate of change of current, dl/dt (called sLew rate) to 
V s /L (V s = DC supply in volts). The delay problem is more severe in systems with 
PWM amplifiers. An experience was recently reported and well explained by Brad- 
field, et al. (1987). To take into account this effect of inductive phase lag on 
the system dynamics, an approximate linear transfer function is 

G = I/E = l/[ ( L/R )S + 1] (14) 

In the above, we neglected the component of back emf, due to variations of in- 
ductance as a result of changes in the air gap. For clarity, the maximum amplitude 
of the transfer function G has been assumed to be one. According to equation 14, it 
is only close to one at U) < R/L. Above R/L, the control current is attenuated and 
delayed. The phase advance by the controller is cancelled by the phase lag in the 
amplifier system. This means negative AMB damping can occur in the high frequency 
range, and the control of flexible rotor criticals can be in jeopardy. 

ROTOR-AMB SYSTEM DYNAMICS 

AMB dynamic stiffnesses are generally lower than ball bearings or oil-film 
bearings. On a typical critical speed map, as shown in Figure 7, one can expect 
that the first two criticals will have rigid mode shapes, and they can be easily 
controlled. The higher modes (particularly the third and the fourth modes) with 
bending mode shapes must be given careful consideration in the design of the AMB 
controls. The first requirement is to assure the existence of positive damping for 
these modes, if the AMBs are responsive to the modal frequencies. Although oil-film 
bearings always provide positive damping, there is no guarantee of this for AMBs. 
Because of the inductance mentioned above, the control current at the high frequen- 
cies may lag behind the displacement measurement. If the operating speed is well 
below the bending critical, one can reduce the lagging current to a minimum at that 
frequency through electronic means. Then, the AMB would not be responsive to the 
modal resonance. 

If the operating speed is above the bending critical, it is imperative to pro- 
vide enough leading current at the appropriate frequency. Even once the proper con- 
trol response has been achieved, the ability of the AMBs to control rotordynamic 
response will be significantly affected by rotor mechanical design. Care must be 
taken in mechanical design to ensure that the AMB locations are not at nodes of the 
vibrational modes of interest. 

This system design challenge can be aided by performing conventional 
rotor-bearing dynamic analyses, such as calculating undamped critical speeds, un- 
balance response, and stability, which have been performed by MTI as part of the 
magnetic bearing design analysis. Hustak, et al. (1985) have performed similar ana- 
lyses for two compressor rotors supported by AMBs. They have shown measured AMB 
stiffness and damping coefficients as functions of excitation frequency. A modified 
but more rigorous rotordynamics approach, which has the inherent advantage of inte- 
gration of controller dynamics into the system equations, is described below. 

The state vector of the conventional rotor model is extended to include the 
state variables of the AMBs. The controller dynamics of each AMB axis, such as 
shown in Figure 5, and the amplifier dynamics (equation 14) can be represented to be 
a set of first-order linear differential equations in terms of the extended states. 
The coupling terms between the rotor model and the AMB model exist in the AMB dynam- 
ic force, which is represented by equation 8. The resultant electomechani cal model 
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using finite element formuLation for the rotor part, can be used for eigenva- 
lue/eigenvector and force response analyses. The rotor model input is the same as 
conventional model including sections of shaft with spec! lied inner diameter, outer 
diameter, length, and concentrated mass and inertias. The input for each bearing 
would be the bearing station number, the measurement station number, and the parame- 
ters Ki, K m , a, b, C d , C v , C e , and L/R. For properly sized AMBs, the feedback 

gains, C d and C v , are the key parameters that can be optimally determined by per- 
forming the eigen and response analyses. This approach to integrated rotordynamic 
analysis has been applied to typical magnetic bearings developed at MTI . 

CONCLUSIONS 

The essence of control and the properties of active magnetic bearings have been 
described and quantified in a language that is more adaptable to mechanical engi- 
neers, In doing so, the following viewpoints on AMB development have been empha- 
sized : 

1. AMB should be treated as a locally controlled and independent device simi- 
lar to other types of bearings. 

2. AMB control axes should be made independent from each other, with identical 
but tunable electronics. 

3. Parallel circuits for stiffness and damping control should be implemented 
to permit more versatility in the application of these bearings. 

4. Pitfalls exist in controlling the rotor bending critical modes; the reason- 
ing, and design guidelines can be explained easily in mechanical terminol- 
ogy- 

5. Modified rotordynamics analysis method with a state vector extended to in- 
clude the AMB state variables are needed for rigorous system performance 
predict ion. 

These mechanical engineering viewpoints are long overdue in gaining prominence 
and influencing the development of AMB technology. 
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NOMENCLATURE 

Active magnetic bearings 
Pole cross-sectional area 
Phase-lead zero parameter 
AMB damping coefficient 
Phase-lead pole parameter 
Radial air gap 
Proportional feedback gain 

Derivative control gain; Phase-lead feedback gain 

Integral feedback gain 

AMB controller output voltage signal 

AMB perturbed magnetic force 

Power amplifier gain 

Self-inductance of electromagnets 

AMB axial length 

Bias current 

Control current 
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/-T" 

AMB stiffness coefficient 
Current stiffness 

Magnetic stiffness due to bias currents 

Dynamic mass 

Number of winding turns 

Integrator output 

Resistance in current path 

Laplace variable 

Time 


AMB static load 

AMB journal displacement in Y-direction 
Phase-lead circuit output 
Excitation frequency 
AMB critical frequency 
integrator cut-off frequency 

d 2 /dt 2 
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Fig. 3 Elimination of Low-Frequency Drift through Integral Control 
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